A phenomenological model of the superconducting state of the Bechgaard salts 
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We present a group theoretical analysis of the superconducting state of the Bechgaard salts, e.g., 
(TMTSF)2PF6 or (TMTSF)2C104. We show that there are eight symmetry distinct supercon- 
ducting states. Of these only the (fully gapped, even frequency, p-wave, triplet) 'polar state' is 
consistent with the full range of the experiments on the Bechgaard salts. The gap of the polar state 
is d(k) oc {tpuk, 0, 0), where tpuk may be any odd parity function that is translationally invariant. 
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One of central challenges facing theoretical physics is 
the full microscopic understanding of unconventional su- 
perconductivity. An important first step towards this 
daunting task is the identification of the correct phe- 
nomenological description of the relevant materials. In- 
deed, our current understanding of the cuprate , heavy 
fermion 0, ruthenate 0, colbaltate fl, and quasi- two- 
dimensional organic superconductors depends on phe- 
nomenological descriptions as in each of these cases there 
is no widely agreed upon microscopic description of the 
superconductivity. However, despite long standing ev- 
idence 0, 01 of unconventional superconductivity in the 
Bechgaard salts and theoretical proposals of triplet states 
[^, the correct phenomenological description of the su- 
perconducting state has not, until now, been identified. 

In this Letter, we perform a group theoretical clas- 
sification of all of the possible superconducting states 
in the Bechgaard salts that respect translational sym- 
metry. This shows that there are only eight symme- 
try distinct states. By considering the properties of 
these states we show that only one of them is consis- 
tent with the full raiige of thermodynamic measurements 
miSIIlIIlIIlIIlliSlillil that have been performed 
on both (TMTSF)2PF6 and (TMTSF)2C104. This state 
is somewhat analogous to the polar state, first discussed 
in the context of ■^Hc. 

Symmetry analysis: The Bechgaard salts form triclinic 
crystals whose symmetry is represented by the Ci point 
group. Ci contains only two elements, the identity and 
inversion. Thus the point group only differentiates be- 
tween even and odd parity states (which we henceforth 
refer to as s-wave and p-wave states respectively) . Note 
that symmetry does not distinguish, say 'd-wave' states 
from s-wave states or '/-wave' states from p-wave states 
as the crystal has neither rotational nor mirror symme- 
tries. (Hence the terms 'd-wave' and '/-wave' are rather 
meaningless in the context of the Bechgaard salts.) All 
the superconducting states unambiguously identified in 
nature thus far, are even under frequency reversal. How- 
ever, this is not required a priori in the superconduct- 
ing state and so we must distinguish between even- or 
odd-frequency pairing the wavefunction of a 

fermionic system must be antisymmetric under the ex- 



change of all labels, the allowed states are then: even- 
frequency, s-wave singlet; odd-frequency, p-wave, sin- 
glet; even-frequency, p-wave triplet; and odd-frequency, 
s-wave triplet. 

The gap function of the singlet phases may be writ- 
ten as A(k) = Ty^/^k, where ?/;k may be any function with 
the appropriate parity that satisfies translation invari- 
ance, and rj is the complex Ginzburg-Landau (GL) order 
parameter. Thus there are only two symmetry distinct 
singlet states: the conventional s-wave, even-frequency 
singlet and a p-wave, odd-frequency singlet. 

To describe triplet superconductivity one must intro- 
duce a complex vector gap function, d(k) [l9| . The in- 
terpretation of d(k) is straightforward (at least for uni- 
tary vectors [lOl) as it points along the Sz = projec- 
tion, i.e., perpendicular to the spin of the Cooper pairs. 
For triplet superconductors we must distinguish between 
weak and strong spin-orbit coupling (SOC). If the SOC is 
sufficiently weak we may rotate the spin and spatial de- 
grees of freedom independently; therefore the symmetry 
group of the normal state is ^ = 5*0(3) (8> G (g) U{1) ^ T, 
where 50(3) is spin rotation symmetry, G is the point 
group of the crystal {Ci for the Bechgaard salts), U{1) is 
gauge symmetry and T is time reversal symmetry. How- 
ever, for strong SOC the independent rotation of the 
spin and spatial degrees of freedom is not a symmetry 
of the system. Therefore, the symmetry group becomes 
g = C^-^'^ (g) U{1) (g) T, where the group C^-''^ is identical 
to the usual point group of the crystal except that the 
operations of the point group simultaneously act on both 
the spin and the spatial degrees of freedom. 

In general when SOC is strong we expect that d(k) cx 
(aXk, b'lk, c*Zk), where Xk, Ik, and Zk are arbi- 
trary functions which respect translational symmetry and 
transform like fc^;, ky, and kz respectively under the sym- 
metry operations of the point group However, as 
Ci contains only the identity and inversion the only re- 
striction on Xk, ik, and is that they have the appro- 
priate parity. Thus, d(k) = T]{a^p^,h'^p^,c*ip^), where 
and may be any functions which have the 
required parity and satisfy translation symmetry, and 
the GL order parameter is a single complex number (rj). 
Therefore the GL free energy is Fg — Fn = a|?7p + /?|?7|*. 



Clearly there is only one solution and therefore there is 
only one symmetry distinct spin part of the wavefunc- 
tion for triplet superconductivity when SOC is strong. 
These states are analogous to the BW phase, which is re- 
alised at ambient pressure in "^He, but may be either an 
even-frequency p-wave or an odd-frequency s-wave state. 

Because Ci only contains one-dimensional irreducible 
representations d(k) = rf ■ *k for weak SOC, where fj 
is the complex vector GL order parameter and '^\^ ~ 
'0kj V^)- The GL expression for the free energy is 
thus U.mFs-Fr, = a|77|2 + /3i|,y|4 + /32 Irf • The 
ground state, up to arbitrary rotations in spin-space, is 
77 oc (l,z,0) for /32 > and ?7 oc (1,0,0) for /?2 < 0. 
We refer to these states as the (3 and polar phases re- 
spectively by analogy with '^He The /? phase cor- 
responds to pairing in a single spin channel, while that 
polar phase corresponds to pairing in both equal spin 
pairing (ESP) channels. Note that a representation must 
be at least two-dimensional for the ABM phase (which 
is realised under pressure in '^He and probably in 
Sr2Ru04 0) to be possible and so this phase can be 
immediately ruled out of our consideration of the Bech- 
gaard salts. Both the (3 and polar states may exist as 
either p-wave, even-frequency, triplet states or as s-wave, 
odd-frequency, triplet states. Thus there arc four possi- 
ble triplet states if SOC is weak. 

Properties of the superconducting states: We list all 
eight symmetry distinct superconducting states for the 
group Ci, which represents the symmetry of crystals of 
the Bechgaard salts, in table U Our task is now to de- 
termine the properties of the of these states and to com- 
pare these properties with those found experimentally in 
(TMTSF)2PF6 and (TMTSF)2C104. 

None of the four even-frequency states are required by 
symmetry to have nodes in the order parameter. This 
extremely unusual property for an unconventional su- 
perconductor is not, as is often stated, because of the 
quasi-one-dimensionality of the Fermi surface, but is a 
direct consequence of the extremely low symmetry of the 
Bechgaard salts. Recall that the basis functions may 
be any function with the appropriate parity. Therefore, 
s-wave states have no symmetry required nodes and p- 
wave states are required to vanish only at the origin (F- 
point), which, by symmetry, the Fermi surface may not 
cross. In contrast, odd-frequency pairing states are in- 
trinsically gaple ss |l8l . Specific heat measurements on 
(TMTSF)2PF6 |l4 and thermal conductivity measure- 
ment on (TMTSF)2C104 ^ both show an exponentially 
activated behaviour, suggestive of a nodeless gap. As 
these experiments see a full gap they are inconsistent 
with odd-frequency pairing. However, the NMR relax- 
ation rate, l/Ti, has a power law temperature depen- 
dence 0, 0| . If this power law were assumed to arise 
from quasiparticles then it would be suggestive of nodes 
in the gap. Rostunov et al. have recently shown that, 
in a triplet superconductor, collective spin-wave excita- 



tions can also lead to a power law dependence of 1/Ti. 
This theory may also resolve the puzzle of why the power 
law dependence of 1 /Ti is seen even at temperatures very 
close to the critical temperature 0, ^|, which is not 
expected from nodal quasiparticles. Thus these experi- 
ments may suggest a triplet pairing state. 

A extremely small peak is seen in 1/Ti, just below 
Tc 0, 01 . However, this peak more than an order of 
magnitude smaller than the Hebel-Slichter expected for 
an even- frequency, s-wave, singlet order parameter |22l| . 
This strongly suggests that the even-frequency, s-wave, 
singlet order parameter is not realised in the Bechgaard 
salts. Further evidence for this conclusion comes from 
the observed strong suppression of the superconducting 
critical temperature by disorder The only state for 

which this suppression of Tc by disorder is not expected is 
the even-frequency, s-wave, singlet order parameter |23l |. 

Evidence for triplet pairing comes from the obser- 
vation that the upper critical fields in the conducting 
planes of (TMTSF)2PF6 and (TMTSF)2C104 exceed the 
weak coupling Clogston-Chandrasekhar (or Pauli) limit 
by more than a factor of four [l^ lisl ll^ . The Pauli 
limit occurs when the Zeeman energy penalty for form- 
ing Sz — Cooper pairs exceeds the condensation energy 
gained by entering the superconducting state and applies 
to singlet states and pairs in the Sz = projection of a 
triplet state When SOC is weak the spin part of the 
order parameter is not 'pinned' to the lattice. Therefore, 
the superconductor may minimise its energy by aligning 
d(k) ± H m nil. Thus, the triplet phases for weak 
SOC will always be ESP phases in the reference frame of 
the magnetic field and arc therefore not Pauli limited. 

In contrast the triplet phases for strong SOC are 
'pinned' to the lattice as the symmetry group does not 
allow the independent rotation of the spin and spatial de- 
grees of freedom. When a field which exceeds that Pauli 
limit is applied to system, it will completely suppress the 
pairing in the Sz = channel, i.e., d(k) goes to zero 
in the direction parallel to the field. However, ESP is 
not suppressed, i.e., d(k) remains finite perpendicular to 
the field. Both possible triplet states for strong SOC in 
the Bechgaard salts have finite components perpendicu- 
lar to the conducting plane and thus we do not expect 
them to be Pauli limited. Hence, we do not expect any of 
the symmetry distinct triplet states to be Pauli limited. 
Therefore, while the large critical field is strong evidence 
against singlet pairing it does not differentiate among 
the six candidate triplet states. It is also worth noting 
that calculations suggest that the observed critical field 
is too large to be accounted for by FFLO singlet states 
which break translational symmetry . Such states are 
therefore not considered in detail in this Letter. 

If a superconducting order parameter breaks TRS then 
spontaneous supercurrents will flow around impurities 
and near grain boundaries |2^ . The most sensitive probe 
of these tiny currents is the zero field muon spin re- 
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TABLE I: Summary of the thermodynamic properties of the eight symmetry distinct states allowed for superconductors with 
Ci point groups and comparison with experiments on the Bechgaard salts. We see that only the polar p-wave triplet state is 
compatible with experiment. Xa{0)/Xn is the ratio of the spin susceptibility in the limit T — » to that in the normal state 
above Tc- The basis functions i/i^j^ (i/^gk) may be any odd (even) parity functions which satisfy translation symmetry. The 
symbol '?' in the experimental sections indicates that an experiment has not been performed. 



laxation (ZF-/iSR) rate. Small fields consistent with 
broken TRS (BTRS) have been observed in UPta il, 
Ui_^Th^Bci3 [H, PrOs4Sbi2 113, and Sr2Ru04 El. 
Importantly, as the currents due to a superconducting 
order parameter with BTRS are extremely small the 
magnetic fields they generate can be suppressed by very 
small longitudinal fields (50 G is sufficient in Sr2Ru04 
[sif'). Luke et al. 0| measured the ZF-^SR rate in 
(TMTSF)2C104 and did not find any increase in the 
relaxation rate to within their experimental resolution 
(« 25 G). This indicates that the superconducting state 
of (TMTSF)2C104 does not break TRS. Both the even- 
and odd-frequency pairing /3 phases are inconsistent with 
this experiment. 

The spin susceptibility, y.^jT ) st rongly distinguishes 
between different triplet states [l9|. When an S'^ = 
pair forms it no longer contributes to the spin suscepti- 
bility and thus for a singlet superconductor or a triplet 
superconductor with d(k)||H, Xs{T) ^ as T ^ 0. On 
the other hand ESP does not affect the susceptibility. 
Thus, Xs{T) does not change upon passing thou gh T r for 
ESP-only states such as the ABM phase of ^He For 
triplet states which contain both ESP and 5*2 = pairs 
the decrease in Xs{T) is proportional to the fraction of 
Sz — Q pairs. As discussed above, the symmetry distinct 
states for weak SOC are all ESP states and thus we ex- 
pect Xs{T)/xn = 1 for all T < where Xn is the spin 
susceptibility of the normal state. In contrast, in both 
the possible states for strong SOC contain pairing in all 
three Sz projections. Thus, 



Xn 



1 - 



(i*ki^ 
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where, {. . ■)fs indicates the average over the Fermi sur- 
face, i G {a, 6, c}, and the superscript on the suscepti- 
bility indicates the orientation of the field 0|. If the 



averages over the three basis functions arc the same, as 
they arc in the BW phase of "^He, Xs (0) /xn = 2/3 for all 
field orientations. No decrease in x(r) is detected below 
Tc in (^TSF)2PF6 with the field aligned along either 
the a 113 or b' [H axes. x^(0)/Xn = Xs(0)/Xn = 1 if 
and only if and vanish everywhere on the Fermi 
surface. As we do not, in general, expect ip^ and to 
vanish everywhere on the Fermi surface the BW state is 
incompatible with the measured Knight shift. 

We summarise the properties of the eight symmetry 
distinct superconducting phases in tabic It can readily 
be seen that the only state consistent with all of the ex- 
periments is the 'polar' p-wave triplet state realised for 
weak SOC which is specified by d(k) oc (^/'k,0,0). It 
is worth stressing the very small number of assumptions 
that have been made to reach this conclusion: (i) trans- 
lational symmetry is not violated by superconducting 
state of the Bechgaard salts; (ii) there is not an acciden- 
tal vanishing of two independent basis functions if SOC 
is strong; and (iii) that the superconducting states of 
(TMTSF)2PF6 and (TMTSF)2C104 have the same sym- 
metry, i.e., that the materials are related by 'chemical 
pressure', (i) is supported by calculations of the criti- 
cal field for the FFLO state in these materials |23 , but a 
triplet analogue of FFLO cannot formally be ruled out at 
this stage. However, it is difficult to see how such a phase 
would be stabilised, (ii) can be tested experimentally: 
Xs(r) ^ as T ^ with H||c* for the BW state with 



(iV'kP)fS = 0; whereas for the the polar 



state Xs{T) = Xn for H||c*. (iii) is easily tested experi- 
mentally; in particular measurements of the Knight shift 
in (TMTSF)2C104 and the ZF-/^SR rate in (TMTSF)2- 
PFg would complete the set of measurements required 
to uniquely determine the superconducting state of each 
material individually. For sufficiently small magnetic 
fields SOC will be strong therefore our analysis predicts 
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that a Fredericks transition 0, [s^ from the BW state 
to the polar state occurs at extremely low fields. The 
different states correspond to different broken symme- 
tries therefore measurements of collective modes [l^l |21| 
could also provide confirmation of our identification of 
the superconducting state. 

Previous theoretical work: Early theoretical work only 
focused on whether the superconductivity was singlet or 
triplet and did not propose a specific triplet state 0- 
The state discussed by Lebed and coworkers [l^ as- 
sumes strong SOC and is therefore a special case of the 
BW phase. Lebed et a/.'s state has (|'0kP)-P''S — but 
{\^k\'^)FS 7^ 0. Eq. (P) shows that this theory will pre- 
dict a large decrease in the spin susceptibility below Tc 
when H|ja. This prediction is clearly contradicted by ex- 
periment [loj and therefore this theory can be ruled out. 
Duncan et al. |33l | have discussed the symmetry distinct 
triplet states in an orthorhombic {D2h) crystal. As the 
Bechgaard salts are triclinic and the angles involved are 
rather large this is not a good approximation. Neverthe- 
less the state they propose [a ^px' state, d(k) oc (/cj;, 0, 0) 
and weak SOC] is a special case of the polar state we have 
shown to be the actual superconducting state. However, 
the ^Px'' state has an accidental node in the plane kx 
As this plane does not cut the Fermi surface the physi- 
cal properties of the ^px ' state are rather similar to polar 
state. However, we stress that the node in the 'p^;' state 
is not required by symmetry and will raise the energy 
of the state, therefore this node is unphysical. Shima- 
hara j^^l proposed that a singlet state is found at low 
field and a triplet state is found in large fields. As inver- 
sion is a symmetry of the crystal, such a change must be 
accompanied by a phase transition. This has not been 
observed and so this theory does not seem compatible 
with experiment. 

In conclusion we have shown that their are eight sym- 
metry distinct superconducting states in monoclinic crys- 
tals with the Ci point group. Of these only the p- 
wave, even-frequency, polar state [d(k) oc (i/'Mk, 0, 0)] is 
consistent with the full range of experiments on both 
(TMTSF)2C104 and (TMTSF)2PF6. There is not yet 
sufficient experimental evidence to determine the super- 
conducting state of cither these materials individually, 
but the chemical pressure hypothesis suggests that the 
polar state is also realised in these materials. Finally, 
we note that the same symmetry analysis applies to the 
Fabre salts. 

This work was motivated by conversations with A. Ar- 
davan. I thank J. Annett, S. Blundell, J. Fjffirestad, 
R. McKenzie, F. Pratt, and J. Varghese for enlighten- 
ing conversations. This work was funded by the ARC. 
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